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ABSTRACT 


The  high  conductivity  of  the  ionosphere  in  large  measure  deter- 
mines the  propagation  of  terrestrial  radio  waves.      The  ionosphere  can 
be  treated  theoretically  as  a  magneto-ionic  plasma  composed  of  elec- 
trons,   ions,  and  neutral  particles.      Certain  macroscopic  properties  of 
these  particles,    such  as  the  complex  index  of  refraction,    can  be  de- 
duced from  the  microscopic  particle  statistics  of  ionized  gases. 
These  properties  can  then  be  applied  directly  to  Maxwell's  equations 
to  determine  the  reflection  and  transmission  process,    provided  a 
suitable  composition  for  the  model  plasma  has  been  found. 

The  composition  of  the  real  ionosphere  is  not  uniform,    especial- 
ly in  the  vertical  direction.    Such  nonuniformity  can  be  treated  theoreti- 
cally with  tne  aid  of  a  flexible  model  ionosphere  employing  the  notion 
of  a  continuously  stratified  plasma. 

An  extensive  survey  of  both  the  index  of  refraction  of  typical 
electron-ion  plasmas  as  models  of  the  ionosphere  together  with  the 
reflection  coefficients  is  made  in  this  paper  to  determine  the  signifi- 
cance of  the  notion  of  heavy-particle  electron- collision  frequencies 
proportional  to  energy. 

The  introduction  of  electron-ion  collision  frequencies  with  a 
linear  dependence  upon  the  average  value  of  the  Maxwellian  particle 
velocity  distribution  results  in  some  interesting  changes  in  the  detail 
of  the  reflection  and  transmission  coefficients  at  low  frequencies.     It 
can  be  noted,   however,    that  the  reflection  coefficients  determined  by 
application  of  the  continuously  stratified  notion  to  particular  models 
of  the  ionosphere  are  not  drastically  changed  by  the  introduction  of 
such  electron  collisions. 


IV 


ON  PLASMA  COLLISION  FREQUENCIES  PROPORTIONAL  TO  ENERGY 
IN  THE  RADIO  WAVE  REFLECTION  AND  TRANSMISSION  PROCESS 

by 

J.    Ralph  Johler 

and 
John  D.    Harper,    Jr. 

1.     INTRODUCTION 
Both  the  reflection  and    transmission   process    of   a   radio 
wave  at  an  electron-ion  plasma  have  considerable  theoretical  and 
practical  importance.     Thus,    theoretical  treatment  of  the    reflec- 
tion process  is  necessary  to  describe  the  field  of  a  radio  wave  trans- 
mitted around  the  earth  by  reflections  from  the  ionosphere.     Also, 
the  field  of  a  radio  wave  at  a  space  vehicle  either  in  or  above  the 
ionosphere  is  a  timely  consideration  for  the  theoretician. 

The  treatment  of  a  continuously  varying  medium  (varying  in  one 
direction,    such  as  the  vertical  z-direction,   Fig,    1)  by  approximating 
the  medium  with  one  or  more  slabs  of  uniform  composition  has  been 
utilized  by  many  authors,    such  as  Hines  [1951],   Wait  [1953],    Ferraro 
and  Gibbons  [1959],   Brekhovskikh  [I960],    Johler  and  Harper  [1962], 
Johler  [  1962]  .     This  paper  applies  such  techniques  to  a  magneto-ionic 
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plasma  carrying  the  process  to  the  limit  until  the  increase  in  the 
number  of  slabs  and  the  decrease  in  the  thickness  of  such  slabs  result 
in  a  convergence  to  a  stable  reflection  or  transmission  coefficient. 
Recently  the  work  of  Phelps  [  I960]    has  developed  the  notion  of 

electron- ion  collision  frequencies  with  a  linear  dependence  on  the 
energy  exponent  of  the  Maxwellian  particle-velocity  distribution. 
This  has  raised  grave  questions  as  to  the  validity  of  the  use  of  a  con- 
stant collision  frequency  in  the  classical  magneto- ionic  theory,    espe- 
cially at  low  frequencies.     This  paper  introduces  the  monoenergetic 
electron  collisions  (electron  collisions  proportional  to  energy)  into 
the  analysis  of  nonuniform  plasma  reflections  and  transmissions  so 
that  a  comparison  can  be  made  with  the  classical  magneto-ionic 
theory,    especially  at  low  frequencies.     An  extensive  survey  of  the 
index  of  refraction  of  an  electron-heavy  particle  magneto-plasma  is 
made  together  with  appropriate  reflection  coefficients  to  determine 
the  significance  of  electron-ion  collision  frequencies  proportional  to 
energy. 

2.     THE  COMPLEX  INDEX  OF  REFRACTION  OF  THE  PLASMA 
The  complex  index  of  refraction  of  a  non-uniform  medium 
varying  in  one  direction,    such  as  the  vertical  z-direction,    Fig.    2, 
will  be  examined  in  some  detail  in  this  paper.     This  is  accomplished 
by  subdividing  the  medium  into  slabs  of  uniform  composition.     It  is 
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therefore  necessary  initially  to  determine  the  complex  index  of 

refraction  of  each  such  slab. 

The  index  of  refraction  of  a  uniform  magneto-ionic  plasma  can 

be  examined  in  detail  with  the  aid  of  the  plane  coordinate  system 
illustrated,    Fig„   1„     The  xy-^plane  is  situated  at  the  lower  boundary 
of  the  model  plasma  defined  as  a  position  in  space  at  which  the  electron 
density  is  nil    (N  =  0)  and  above  which(z  >  0)  the  ionization  varies  as 
some  arbitrary  function  (or  profile),    N(z),    N(h),   where    h  in  practice 
is  the  height  above  the  terrestrial  sphere  in  radio  wave  propagation 
problems.     Also,    at  all  z  <  0,    the  ionization  is  assumed  to  be  negligible. 

The  mathematical  description'of  the  propagation  in  a  magneto- 
ionic  medium  can  be  formulated  with  the  aid  of  a  plane  wave,     E., 

i 

incident  from  the  region  below  the  xy-plane  (z<  0,    Fig.   1)  at  an  angle 
of  incidence  with  the  vertical  z-direction,    <j>., 

E.  =  |E.  I   exp  [i  (cot  -  -    y]  D)]  (1) 

1       '     l '  c 

where  the  field  is  assumed  to  be  a  time-harmonic  wave  oscillating  in 

o 

time,    t,    at  a  frequency,     f=  co/ 2tt,   with    c    the  speed  of  light,  c  ~  3(10   ) 

meter/ second,    and    r\    the  index  of  refraction,   which  for  a  vacuum; 

Ti  =  r\     =1,    and  D  =  ti     D  =  x  sin  d>.   sin  <j>      +  y  sin  <t>.  cos  <t>     +   z  cos  <t>..  '    ' 
o  'o  i  a       '  l  a  i 

The  magnetic  azimuth,     $    ,     is  reckoned  clockwise  from  magnetic 

ct 

north,  Fig.  1,    i.e.,   from  the  yz-plane,   which  contains  the  static 


t  Rationalized  MKS  units  are  used  throughout  this  paper. 
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magnetic  field  vector,     W ,   oriented  at  a  dip  of  inclination  angle 
relative  to  the  horizontal,   I.     A  resultant  wave,     E  ,    transmitted  into 
the  model  ionosphere  magneto-ionic  plasma  (z  >  0)  is  then  assumed 
to  have  the  form, 

Et=    |Ej    exp[i(cot  -^tj'D)],  (3) 

where  in  the  plasma, 

T)  D  =  x  sin  <t> .  sin  <t>     +  y  sin  $ .  cos  <t>     +   z  £,,  (4) 

i  a  l  a 

in  which    t,    is  in  general  a  complex  number,    the  value  of  which  will 

depend  upon  the  altitude,    z,    and  on  £,  =  C,(z),    since  the  complex  index 

of  refraction,  t|, 

2         2  2 

r)     =  C     +   sin    $.,  (5) 

varies  with  the  electron-ion  composition  which  in  turn  changes  with 
altitude. 

The  complex  index  of  refraction,  T),  can  be  determined  locally  at  a 
particular  altitude  by  a  simultaneous  solution  of  Maxwell's  equations, 

VX     E+    HQff  =0, 

VX    H    -  J  -  €    |p-=  0,  (6) 

o  at 

where  E  and  H  are  the  electric  (volts /meter)  and  magnetic  (ampere- 
turns/meter)  fields,  J  =  -NeV"   the  vector  convection  current  (amperes/ 
square  meter), 
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and  the  Langevin  equation  of  motion  of  an  electron, 

dV  —  —       -. 

m—  +  mg  V     +   |i  e  (V  X  W)+    eE  =  0,  (7) 

in  which  m  is  the  electronic  mass,    e    the  electronic  charge,    p.       the 

permeability  of  space  (jj     =  4n(10      )  henry/meter),    &  the  vector  static 

magnetic  field  (ampere-turns/meter),    V  is  the  vector  velocity,    and  g 
is  the  collision  frequency  parameter  (seconds        ).     The  elimination  of 
the  vectors  J  and  H  from  the  simultaneous  solution  of  equations  (6,    7) 
results  in  the  matrix  equation,   using  the  Cartesian  coordinates  of  Fig.  1, 


ail        al2         al3 


a21       a22        a23 


a31       a32        a33 


x 


=    0, 


(8) 


where   the  coefficients  will  be  defined  in  Appendix  I. 

If  the  field  E,    equation  (8),    exists  the  determinant  of  the  matrix  must 
vanish.     Expansion  of  the  determinant  yields  a  quartic  equation  in 


£,     [Booker,    1939],    [Johler  and  Harper,   1962], 
a4l4+   a3C3  +   a2£2+   a^  +   aQ    =  0, 


(9) 


with  coefficients  also  defined  in   Appendix  I. 


The  Booker  [1939]    quartic  equation  (9)  is  quite  significant, 

since  the  four  complex  roots  of  such  an  equation  describe  four  indexes 

i   r 
of  refraction,    r\  =  r\   '  '      corresponding  to  upgoing  (i)  with  negative 

imaginary  part  (Im  t\  <  0)  and  downgoing  (r)  with  positive  imaginary 
part  (Im  r\  >  0)  together  with  ordinary  (o)  and  extraordinary  (e)  for 
each  upgoing  or  downgoing  magneto- ionic  propagation  component.    The 
intensity  and  phase  of  these  propagation  components  are  dependent 
upon  direction  of  propagation  relative  to  the  direction  of  the  static- 
magnetic  field,    &,    (Fig.  1)  which  acts  on  the  convection  currents,    J, 

as  the  Lorentz  force,    equation  (7), 

|i 

-^(Jxl)  =  ne  (VX    U). 
IN  O 

Such  a  medium  of  propagation  is  said  to  be  anisotropic.     Indeed,    the 

ionosphere  which  surrounds  the  earth  is  one  such  plasma  medium. 

3.     BOUNDARY  CONDITIONS 

The  detailed  structure  of  the  flexible  plasma  model  is  illustrated, 

Fig.   2,    as  a  stack  of  plasma  slabs  of  arbitrary  thickness    (except 

for  the  topmost  slab  of  thickness  z     =  °°).   The  number  of  such  slabs, 

P 

p,     and  the   thickness    is   quite   flexible,     since   the   idea   of   this    analy- 
sis   is    that   the   measured    electron  density- altitude    and    collision 
frequency- altitude    (N(z)    and    v(z),     respectively)    profiles    can  be 


approximated  to  any  desired  accuracy  by  decreasing  z     and  increasing 
p  simultaneously  until  a  stable  reflection  process  is  obtained. 

A  constant  electron  density,    collision  frequency,    and  static  mag- 
netic field  with  respect  to  altitude,    z,    is  of  course  assumed  for  each 

slab,    z    ,    and  associated  with  each  such  slab  a  set  of  four  roots,    Q  =  Q    , 
n  n 

is  found  to  exist.     Two  of  the  roots  will  exhibit  a  negative  imaginary 
part  (Im  C,  negative),    corresponding  to  an  upgoing  propagation  compo- 
nent (+  z  direction,    Figs.    1  and  2).     Also,    two  of  these  roots  will    ex- 
hibit a  positive  imaginary  part  (Im  Q  positive),    corresponding  to  a 
downgoing  propagation  component  (-z  direction,    Figs.    1  and  2).     Ex- 
cept for  the  topmost  slab,    it  is  necessary  to  consider  both  upgoing  and 
downgoing  components  in  this  analysis. 

The  treatment  of  a  continuously  varying  medium,    i.  e.  ,    one  which 
varies  in  the  vertical  z-direction,    Figs.    1  and  2,    with  one  or  more 
slabs  of  uniform  composition  was  presented  in  a  previous  paper  [Johler 
and  Harper,    1962]  .     Indeed,    the  technique  was  carried  to  the  limit  at 
which  the  number  of  slabs,    p,    for  each  calculation  depends  upon  the 
computation  precision  required  and  the  particular  values  of  the  electric 
and  geometric   parameters.       Reflection   and    transmission   coefficients 
were  calculated,     and  the    stratification   process    was    continued   for 
model  electron  density   and    collision   frequency   profiles   until   a 
stable    reflection   or   transmission  was    obtained.       It  was  therefore 
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demonstrated,    that  subject  to  the  provision  that  a  convergence  could 
be  demonstrated  numerically,    the  stratified  model  could  be  employed 
to  represent  a  continuum,    hence  the  notion  of  continuous  stratification. 

The  distinction  between  an  ordinary  and  an  extraordinary  propa- 
gation component  is  ambiguous  and  obscure  in  the  nonuniform  iono- 
sphere.    Indeed,    if  only  the  reflection  coefficient  is  of  interest  the 
distinction  is  immaterial.     However,    if  the  progress  of  the  wave   up- 
ward through  the  ionosphere  is    to   be    examined   in  detail,     the  dis- 
tinction should  be  maintained  to  avoid  a  confusing  analysis.     This 
is    accomplished   by  an   examination  of   the   form   of   the   index  of 
refraction  function  with  respect  to  frequency  and  altitude  (or  electron 
density  and  collision  frequency  which  vary  with  respect  to  altitude). 
Thus,    the  complex  index  of  refraction,    equation  (5),    is  detailed  for 

each  frequency  and  slab,    z    ,    r\  =T\   ,    and  the  upgoing  (i)  ordinary  and 

n  n 

extraordinary    r\ .  and  the  downgoing  (r)  ordinary   and    extra- 

1,    o,    e 

ordinary   ri  function  continuity  is  examined  in  detail   as    a  func- 

r,    o,    e 

tion  of  frequency   to  determine   the    crossover   point   of   the   function 
for  each  slab  or    electron  density   under    consideration.       The  Re  r\ 
was    employed   in   this    analysis.       Thus,     above   the    crossover   point 
the  ordinary  wave  was  considered  to  be  the  greater  of  the  two  roots, 
Re  T\  >Re  r\  ,    and  below  the  crossover  point  the  ordinary  wave  was 
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considered  to  be  the  lesser  of  the  two  roots,   Re  ti     <  Re  r\   .     At 

o  e 

precisely  the  crossover  point,    the  two,    of  course,    are  identical, 

-n     =  r\   .     This  point,   for  examples  with  low  electron  density  examined 
o         e 

in  this  analysis,   was  found  to  be  below  the  plasma  frequency, 
f      =  go    /2tt.     The  absolute  distinction  between  the  two  propagation 
components  remains  quite  arbitrary,    but  the  analysis  must  be  con- 
sistent between  each  slab  and  consistent  within  each  slab  for  upgoing 
and  downgoing  propagation  components  to  simplify  the  detailed  study 
of  the  wave  progress  in  the  ionosphere.      The  almost  random  possi- 
bility as  to  order  in  which  the  four  components  emerge  in  the  numer- 
ical solution  [Johler  and  Walters,    I960]  of  the  quartic  equation  (9) 
(the  Muller  [1956]  iterative  procedure  was  employed)  makes  necessary 
this  logical  procedure  for  sorting  out  the  roots  thus  obtained. 

The  boundaries  of  each  slab  are  introduced  as  an  expression 
for  the  continuity  of  the  tangential    E     and    H   fields  and  the  normal  H 
fields  at  each  boundary,   Fig.    2,    of  the  model  plasma.       This    is 
accomplished  by  equating  the  field  immediately  above  and  immediately 
below  each  boundary  which,    after  considerable  ado,    can  be  expressed 
in  concise  form  as  the  matrix  equation  (10), 
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where  the  elements      a. 


can  be  deduced  as  a  consequence  of 
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these  boundary  conditions  with  results  as  shown  in  Appendix  II, 


provided  the  reflection  coefficients,     T,    and  the  transmission 


coefficients,     U,    T  ,   are  defined, 
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where,  Hg.    2,     n  =  1,     2,     3.  .  .  ,     p-1,     p,     and  the  four    T"    are 
defined  as  a  transmission  coefficient    into  a  topmost  slab  of  infinite 

extent,    z,    and  zero  electron  density,     N  =  0,    f]     =  n     =  1. 

p         o 

The  reflection  coefficients,     T      ,    T        ,    T        ,    T  ,    describe 

ee         em        me        mm 

the  reflection  into  the  region  below  the  model  plasma,     z  <  0.     Thus, 

T        refers  to  vertical  electric  polarization  of  the  incident  wave  and 
ee 

a  corresponding  vertical  electric  polarization  of  the  reflected  wave. 

T  refers  to  the  generation  of  the  abnormal  component  (horizontal 

em 

electric  polarization)  by  the  incident    vertical  electric  wave. 

Similarly,     T  refers  to  vertical  magnetic  incident  wave  and 

mm  ° 

vertical  magnetic  reflected  wave,    and    T  refers  to  a  correspond- 

me 

ing  abnormal  component  (horizontal  magnetic). 

The  transmission  coefficients,    U,    refer  to  transmission  at 
the  particular  point  in  the  ionosphere  under  investigation  to  whiGh 
the  wave  has  penetrated.      There  are  four  propagation  components 
and  two  polarizations  and,    hence,    eight  types  of  transmission  coeffi- 
cients:    upgoing  and  downgoing,    ordinary  and  extraordinary,    vertical 
electric  and  vertical  magnetic. 

Thus,    the  nature  of  the  wave  reflected  from  the  plasma,    the 
nature  of  the  wave  transmitted  through  the  plasma,    and  indeed  the 
nature  of  the  wave  progressing  within  the  plasma  are  completely 
described  by  the  coefficients  determined  by  this  analysis. 


-12- 

4.      CLASSICAL  MAGNETO -IONIC  THEORY 
The  Booker  [1939]  quartic  equation  C9)  is  satisfied  by  the 
famous   Appleton-Hartree  formula  [Appleton,   1932],   [Hartree,   1929, 
1931], 

n   =«i /    2S  1  2 <12> 

2(1-—) sm     w   ±       — -  sin     \\>  +   4  — r-  (1  •  — )     cos      oi. 

s  2  \       2  2  s 

s  v    s  s 


2  2  2 

where,     £>     s    "H     -  sin       <j>.,  (13) 

provided, 

t  1 

cos  di  =  -  —    sin  I  +  —    sin  <j>.     cos  <t>     cos  I,  (14) 

r|  7]  i  a  •        1 


or, 


sin  d/   =    — 


2       „  2         2  .  2 

ri     -  £,      sin     I  -  sin  <j>.  cos  <j>      (sin  <t>.  cos  0     cos     I 

1  a  1  a 


-  Zt,  sin  I    cos  I) 


This  formula,    employing  real    v|j,  has  become  identified  with  the 
classical  magneto-ionic  theory-     This  solution  becomes  identical 
with  the  Booker  quartic  with  complex    \\i  [Johler  and  Walters,   I960], 
i.  e.  ,   the  propagation  direction  relative  to  the  earth's  magnetic  field 
is  in  general  complex,     i|i  =  Re  ip  +   i  Irn  i|j. 


(15) 
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5.     ELECTRON  COLLISION  FREQUENCIES  PROPORTIONAL 

TO  ENERGY 

The  particle  statistics  of  the  electron-ion  gas  with  superposed 
electrodynamic  and  magnetostatic  fields  has  been  investigated  in 
some  detail  since     the  classical  magneto -ionic  theory  was  intro- 
duced (1932).     Most  of  these  studies  have  arrived  at  macroscopic  pro- 
perties from  the  well  known  Maxwell-Boltzmann  theory  by  means  of 
approximations  of  the  Boltzmann  equation.     The  Langevin  equation  (7) 
for  constant  collision,     g  ~  v,    is  one  such  approximation.     The  re- 
examination of  the  classical  magneto-ionic  theory  in  recent  years 
can  be  traced  to  the  work  of  Huxley    [1937,   1938,    1940],    Pfister  [1955], 

Jancel  and  Kahan  [1955],    Dingle,   Arndt  and  Roy  [1956],   Molmud  [1959], 
Allis  [1956]  ,     and   Sen  and  Wyller  [I960]  .     Particular  emphasis, 

based  on  experimental  work  with  gases  in  the  laboratory,    has  been 
placed  by  Phelps  [I960]    on  the  dependence  of  the  collision  process 
upon  the  energy,    au        ,    where  u  =  kT/e  (k  =  Boltzmann's  constant, 
T  =  temperature,    degrees  kelvin,    and  e  =  electronic    charge),    with 
special  emphasis  on  the  linear  dependence,   j    =  2,    on  the  collision 
frequency,    V  =  v(u)  =  au,    where  a  is  a  constant.     This  linear  depend- 
ence of  the  collision  frequency,    v(u) ,    upon  energy  has  been  called  a 
monoenergetic  electron  collision  process.     Phelps  [  I960]     anticipated 
grave  consequences  as  to  the  general  validity  of  the  classical  magneto- 
ionic  theory.      This  paper,    therefore,  has  introduced  collision 
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frequencies  proportional  to  energy  into  the  analysis  of  the  reflection 
coefficient,    and  the  consequences  of  such  a  theory  are  examined. 

Phelps  [i960]  has  developed  integrals  of  the  electron  velocity 

or  conductivity  tensor.     The  integration  over  a  Maxwellian  energy 

2 
distribution  expressed  in  terms  of  energy,  eu  =  i  m  V  ,    or,     f     = 

-i. 
[e/ir  k  T]        exp  [-eu/k  T], involves  momentum  transfer  collisions  with 

gas  molecules  proportional  to  the  energy  u,     v    =  v(u)  and  can  be 
accomplished  for  a  continuous  electromagnetic  wave  of  angular  fre- 
quency,   oj  =  2tt  f. 

The  notion  of  a  complex  parameter,    g,    in  the  Langevin  equation 
(7)  was  introduced  by  Molmud  [1959].     Thus,    in  the  classical  theory, 
g  ~  v.   In  the  monoenergetic  theory  the  parameter    g  =  f(v(u))  is  a 
complex  and  frequency-dependent  number.      The  implication  of  the 
classical  magneto-ionic  theory  that  such  a  quantity  can  be  repre- 
sented by  a  real  effective  value  without  a  frequency  dependence  is 
therefore  not  in  general  valid  concept. 

The  derivation  of  the  parameter,    g,    and  its  introduction  into  the 
analysis  is  given  in  Appendix  III. 
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6.     COMPUTATIONS  AND  DISCUSSION 
The  precise  model  for  the  lower  ionosphere  has  been  the  subject 
of  extensive  investigation.     Such  a  model  can  be  presented  as  an  elec- 
tron density-altitude,    N(h)  or  N(z),    Fig.     2,     profile    together   with    a 
collision  frequency- altitude  profile.      Two  such  models,    N(h)  or  N(z) 
profiles,    are  illustrated,    Fig.     3.       The  Houston  Composite  [  Waynick, 
1957]     represents    the    approximate    condition   of   the    ionosphere    at 
daytime-noon.      A  profile  which  was  measured  with  a  rocket  sounding 
by  Sedden  and  Jackson  [  1958]    is  also  illustrated,     Fig.     3.     This  model 
is  representative  of  the  condition  of  the  lower  ionosphere   during   a 
local  auroral-zone  blackout  of  the  type  confined  to  the  auroral  zone 
or  some  segment  thereof.     Measurements  were  not  made  to  establish 
this  model  below  1000  electron/cm3.     However,    various    methods 
have  been  proposed  to  extrapolate  such  models  to  lower  electron  den- 
sities.    The  Chapman  theory  of  ionized -layer  formation  can  be  written, 

Chapman  [  1931,    1939]  ,   Budden  [  1961]  , 

r    z  "  z    -\    n 

(50) 


N  =  N     exp  i 
o       ^2 


z  -  z  z  -  z 

i                   °  o 

1  -  — ^ -  sec  x  exp 


.  1 

where    N     =  N  /(cos  y)2 

o  max  '  A 


z     =    z  -  S  In  (sec  y) 

o  max  A 


S     =  RT/Mg 
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R    is  the  gas  constant,     M   is  the  mean  molecular    weight,     T  is  the 

absolute  temperature,    g  is  the  acceleration  of  gravity,    and  x  is 

the  zenith  angle  of  the  sun.     The  scale  height,  S,    can  in  practice  be 

adjusted  to  fit  measured  profiles  provided    N  and    z  are 

max  max 

determined  by  such  profiles.     The  gaussian  extrapolation  formula  is 

often  employed,    [Budden,   1955],    for  the  daytime-noon  profiles,    in 

which     N     =  N  ,     z     =  z 

o  max         o         max 

N     =  N  exp  [-  (z  -  z  )2A'l, 

max  max     '      J 

where  the  constant    k'    is  adjusted  for  best  fit  of  the  measured  profile. 

This  latter  procedure  was  followed  in  the  ionosonde  blackout  model 

profile  illustrated,    Fig.    3  .       Obviously,    a  variety  of  other  measured 

or  theoretical  profiles  could  be  selected.     These  profiles  were 

employed  primarily  to  demonstrate  the  techniques  described  in  this 

paper.     The  analysis  of  a  variety  of  other  profiles  is  reserved  for 

future  work. 

Numerous  v(z)  profiles,    both  theoretical  and  measured,    based 

upon  the  classical  magneto- ionic  theory,  are  also  available  in  the 

literature.     Some    of   these  illustrated,    Fig.    4,    have  been  deduced 

by  Nicolet  [  1958]  ,    Fejer  [  1955]  ,     Crompton,    Huxley,    and  Sutton 

[  1953]  ,    Gardner  and    Pawsey  [  1953]  ,     and  Sedden  and  Jackson 

[  1958]  .       The    results    of  Sedden   and    Jackson  [  1958]  indicate    the 
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theoretical  values  of  Nicolet  should  be  divided  by  3,   hence  the 

Nicolet/3    v  (z)  profile.      The  Nicolet/3  v(z)  profile  was  selected  for 

the  analysis  presented  in  this  paper.     The  collision  frequency  can  be 

deduced  from  experimental  measurements  at  high  frequencies,    (HF). 

Thus,   for  example,   Kane  [i960]      employed  a  frequency    7.  75  Mc/s 

together  with  the  sixth  harmonic  as  a  reference  and  measured  the 

relative  attenuation  of  the  ordinary  and  extraordinary  waves,    from 

transmitters  in  a  rocket  traveling  upward  through  the  D-region. 

These  attenuations  are  dependent  on    Im   r\         ,    or 

o,  e 


exp 


-  —    Im   Ti        D 
c  o,  e 


The  theoretical  calculations  of  Im  ti  at  7.75  Mc/s  and  2.0  Mc/s 

o,  e  ' 

employing  the    roots   £,  of  the  quartic  equation  (9)  in  equation  (5)  to 

determine   -n  and  using  an  arbitrary  value,    N  =  500,    is  illustrated, 

o,  e 

Figs.    5  and  6,    for  classical  and  monoenergetic  theory.        This 

implies  that  the  numbers  v  for  the  abscissa  of  the  graphs  are   constant 
average  values  of    v     for  the  classical  curve  and  proportional  to 

energy,   u,    or    v  =  v(u)   for  the  monoenergetic  curve.     Upon  selecting 

a  value  on  the  monoenergetic  curves  for  ordinary  and  extraordinary, 

and  keeping    Im   r\    constant  and  moving  over  to  the  classical  curve, 

v  .  4 

one  notes  a  change    — ; — r-    ~    2.  5    for  collision  frequency  values    10     to 

v(u) 

7 
10     approximately,   which  is  in  agreement  with  the  deductions  of 
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Kane  [I960]  .      The  values  of  the  Nicolet/3  curve,    Fig.    4,    based  on  the 

classical  theory  were  therefore  reduced  by  a  factor  of  2.  5  when  used 

in  the  monoenergetic  theory,    V  =  v(u).     It  is  interesting  to  note  that  the 

Nicolet/3  curve  then  agrees  quite  well  with  one  of  the  curves  of  Kane 

[  I960]    deduced  with  the  aid  of  monoenergetic  theory  from  rocket 

measurements. 

Reflection  coefficients,    T       ,    T        ,    T        ,    T  ,    which  can  be 

ee        em        me        mm 

employed  to  determine  the  field  strength  of  a  radio  wave  at  the  surface 

of  the  earth  [  Johler,    1961c]    were  evaluated  with  the  aid  of  the  models, 

Figs.    3,    4,    and  the  continuously  stratified  technique,    Fig.    2. 

The  results  of  the  calculation  of  the  reflection  coefficients  based 

on  the  monoenergetic  theory  are  compared  with  the  results  based  on 

the  classical  theory  as  a  function  of  frequency  for  grazing  incidence 

on  the  ionosphere,    $ .  ~  82°,    (d  ~  1000  statute  miles  for  a  single 

ionosphere  reflection  [  Johler,    1961c]s    using  the  Houston  composite 

model,    Fig.     3.       A   magnetic   inclination   or   dip,     I  =  60    ,    was 

employed,    and    the    results  for  propagation  into  the  east,     magnetic 

azimuth,    <\>     =  90°,    Fig.    7,    and  propagation  into  the  west,    <t>     =  270° 
a  a 

Fig.    8,    are  illustrated.       In  view  of   the    grave    consequences    antici- 
pated by  the  introduction  of  the  electron  collisions  proportional  to 
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energy  it  is  most  remarkable  that  the  curves  of  the  classical  theory- 
differ  only  in  detail  from  those  of  the  monoenergetic  theory.     Indeed, 
for  the  models  employed  in  this  analysis  no  new  gross  phenomenon 
in  the  propagation  were  described  by  the  monoenergetic  theory. 
However,    some  of  the  detailed  changes  are  of  considerable  interest. 
Thus,   for  example,  at  100  kc/s    there  is  a  noteworthy  difference 
between  the  abnormal  components,     T  *»   T        f   based  on  classical 

theory  as  compared  with  monoenergetic  theory.      This  suggests 
comparison  of  loop  and  vertical    antenna   measurements  at  grazing 
incidence  for  first  ionospheric  reflected  wave  from  a  vertical 
radiator  type  transmitter,   for  example,   to  check  the  monoenergetic 

theory  experimentally. 

Similarily,    the  reflection  coefficients,     T,    based  on  monoenergetic 
and  classical  theories  were  compared  with  the  blackout  model  profile 
N(z),    Fig.    3,    employing   the  same  (Nicolet/3)   model   v(h)    profile, 

Fig.    4,    as  shown  in  Figs.    9  and  10  for   propagation  into   the    east 

(cj>     =  90  )   and  propagation  into  the  west    (<t>     =  270  ).    The  geometric- 
a  a 

optical  angle  of  incidence,    <t>.,   was,    of  course,   adjusted  slightly  to 

take  account  of  the  lowering  of  the  ionosphere.     Here  again,   the 

difference  between  monoenergetic  and  classical  are  only  those  of 

small  detail.     The  greatest  effect  seems  again  to  be  in  the  abnormal 

components,     T         =  T 

em         me 
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The  reflection  coefficient  as  a  function  of  angle  of  incidence,    $., 

1 

which  is  related  to   geometric-optical  ray  distance  subtended  along 
surface  of  the  earth  between  transmitter  and  receiver  [Johler,    1961c] 
d,    statute  miles  as   shown  by  the  additional  abscissa,    is  illustrated  for 
the  quiescent  model  using  the  monoenergetic  theory,    Figs.    11,    12. 
The  results  for  propagation  into  the  east,    $     =90°,    Fig.     11,    can  be 

Si 

compared  with  the  results  for  propagation  into  the  west,    <j>     =  270°, 

3, 

Fig.    12.      The  particular  case  of  a  reflection  coefficient  for  vertical 

polarization  of  transmitter  and  receiver,    T        (vertical  electric  inci- 

ee 

dent  wave  and  corresponding  vertical  electric  reflected  wave),    is  illus- 
trated as  a  function  of  distance,    d,    or  angle  of  incidence,    $.,    in 
Fig.   13  as  a  comparison,    employing  the  Houston  Composite   model, 

Fig.    3,     of   east-west,     <t>     =  270      with  west-east,     <t>      =90     propaga- 

a  a 

tion  together  with  a  comparison  of  monoenergetic  with  classical 

3 
magneto-ionic  theories.      Thus,    at  grazing  incidence,    $ .  ^   82°,    d  ~  10 

miles,    I  T       I    monoenergetic  and  classical  results  exhibit  a  greater  field 
1     ee1 

intensity  for  propagation  into  the  east  as   compared  with  propagation 
into  the  west.      This  is  consistent  with  previous  conclusions  based  on 
classical  theory,    only  [Johler  and  Harper,    1962],    [Johler,    1961a], 
[Johler  and  Walters,    I960],    [Johler,    1961c].     For  both  monoenergetic 
theory  and  classical  theory  both  curves  cross  in  the  vicinity  of 
Brewster's  angle  (d  ^   250),    where  at  short  distances  or  near  vertical 
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incidence  the  propagation  into  the  west  is  more  favorable  than  the 
propagation  into  the  east  for  this  particular  orientation  of  the  earth's 
magnetic  field,    I  =  60°,    $     =90°,    270°.      The  fields  predicted  by  the 

cL 

classical  theory  are  less  than  those  predicted  by  the  monoenergetic 
theory  for  propagation  into  the  west.     However,    for  propagation  into 
the  east  there  is  a  cross-over  in  the  curves  as  a  function  of  angle  of 
incidence,    <{>.,    or  distance,    d.     A  distinguishing  feature  of  these  curves 
when  compared  with  ordinary  sharp-boundary-type  Fresnel  reflection 
coefficients  is  the  multiple  Brewster  angles  instead  of  a  single  Brewster 
angle  (between  d  =  90  and  300).      This  can  be  interpreted  physically  as 
multiple  regions  within  the  ionosphere  of  large  contribution  to  the  total 
reflection  coefficient.      Thus,    each  such  region  would  exhibit  a  slightly 
different  Brewster  angle  since  the  electron-ion  composition  of  the 
ionosphere  varies  with  altitude,    z. 

Figures  14-25  illustrate  the  complex  indexes  of  refraction  repre- 
sentative of  low,    medium,    and  high  elect  ron  densities.      The  collision- 
frequency  parameter  also  affects  the  shape  of  the  curves.      Those  chosen 
were  associated  with  the  Houston  Composite  N(h)  profile  in  the  reflec- 
tion coefficient  analyses.      Grazing  incidence,    <t>.  ^  82°  is  assumed, 
and  the  earth's  magnetic  field,    Jt  =  0.  5  gauss,    is  oriented  such  that 
<D     =0°,    and  I  =  60°. 
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In  making  a  comparison  of  the  indexes  derived  from  the 
monoenergetic  theory  with  those  derived  from  the  classical  theory,    perhaps 
the  most  noteworthy  observation  is  the  lack  of  major  differences.     In 
fact,    as  the  electron  density  is  increased  the  two  cannot  be  separated 
graphically.     This  would,    in  large  measure,    account  for  the  only 
minor  changes  in  the  reflection  coefficients  exhibited  by  the  monoenergetic 
theory. 

An  interesting  feature  of  both  the  monoenergetic  and  classical 
index  curves  is  the  tendency,    as  the  electron  density  is  increased, 
for  one  or  both  of  the  extraordinary  propagation  components  to  evanesce  as  a 
function  of  frequency  (i.e.,    the  point  of  evanescence,    Re  r\  =  0)  although 
not  at  the  same  place. 
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7.     CONCLUSIONS 

The  continuously  stratified  concept  together  with  measured  or 
theoretical  N(h),    u(h)  profiles  of  the  ionosphere  can  be  employed  to 
analyze  a  variety  of  conditions  of  the  ionosphere  at  LF.     The  intro- 
duction of  electron  collision  frequencies  proportional  to  energy  into 
the  magneto- ionic  theory  introduces  some  interesting  changes  in  the 
detail  of  the  reflection  coefficients  and  index  of  refraction  of  the 
lower  ionosphere.     However,   the  classical  theory  demonstrates  the 
principal  features  of  the  reflection  process  with  remarkable  accuracy 
in  view  of  the  grave  consequences  anticipated  by  the  introduction  of 
electron  collisions  proportional  to  energy. 

The  introduction  of  electron  collision  frequencies  proportional  to 
energy  provides  a  more  accurate  description  of  the  macroscopic 
properties  of  a  magneto- ionic  piasma.     The  comparative  simplicity 
with  which  such  collision  frequencies  can  be  introduced  into  the  full- 
wave  solution  for  the  ionosphere  warrants  their  use  in  describing  the 
detailed  reflection  or  transmission  process. 


Boulder,   Colorado,   March  29,    1962. 
(Revised  August  8,    1962). 
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APPENDIX   I 


The  elimination  of  the  vectors,    J  and  H,    from  the  simultaneous 
solution  of  Maxwell's  equations  and  the  Langevin  equation  of  motion  of 
the  electron  gives 
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If  the  field,    E,,  exists  the  determinant  of  the  matrix  must  vanish, 
which  determinant  upon  expansion  yields  a  quartic  equation  in  C 
[Booker,    1939]  ,    [  Johler  and  Harper,    1962]  , 

a4£4+   a3C3+   a^2+   a^  +   aQ  =  0, 

where, 
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h      -  -  h  sin  I, 
JL/ 

hT  =  h  cos  I, 
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the  plasma  frequency, 
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and  where    N  =  N       the    median  value  of  electron  density  in  the 


n 


interval   z       corresponding  to  the  uniform  slab  number    n,   Fig.    2 


n 


provided  the  approximation    g  ^  v    or    g     ^  v        for  each  slab  of 
ordered  number     n,    can  be  made,     i.  e.  ,   the  classical  magneto-ionic 
theory  is  assumed  by  a  real  constant  value  of  collision  frequency,    v  . 
Model  electron  density  profiles,    N(z),  and  collision  frequency    profiles, 
v(z),    are  shown  in  Figs.    3  and  4.     Such  profiles  are  considered  to  be 
representative  of  daytime  noon  and  local  auroral  zone  blackout  con- 
ditions of  the  earth's  ionosphere,    and  therefore  have  considerable 
practical  significance  to  radio  wave  propagation  science. 
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APPENDIX  II 


The  coefficients  of  the  matrix  equation  (10)     can  be  written: 
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Field  ratios  in  the  plasma  are  of  considerable  interest. 
Defining, 
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equation  (8)  can  be  written, 
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s    -h     -J 


s    -h 


a^-i 
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2   i2 
s    -h 


w 


here  the  particular  slab    (Fig.   1),    n  =  1,     2,     3,     .  .  .  p-1,    p,    under 


consideration  is  designated  by  the  notation,      C  =  ^         ,P=P       ,    Q  =  Q 
Since  there  are  four  roots  to  the  quartic,    the  field  ratios  equations 

(27,    28)  are     PW,  P(n),  P(n),  P(n),  Q(n),  Q(n),  Q(n)  ,  Q(n),   corres- 
io         le  ro         re         10  le  ro  re 

ponding  to  the  zetas  (£,)  for  upgoing  and  downgoing  (i,  r),    ordinary  and 
extraordinary  (o,  e)  propagation  components  in  the  particular  region, 


z    ,    of  the  electron-ion  magnetic  plasma  under  consideration. 
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APPENDIX  III 


The  velocity,    V,    tensor  can  be  written, 


J  =  _  NeV  = 


Al  AZ  A3 
Bl  B2  B3 
Cl    °2    C3 


E 


y 


(i) 


where 


s(iweQ) 

Al   =       2      2      ' 
s    -h 


Bl   = 


-ihL(i(V 

2      2 

s    -h 


ihT(ia)€Q) 

2  ,  2         ' 
s    -h 


2  2 

ihL(iweQ)  -(h      -  s    )(iwc    )  -h    h    (iwc    ) 

A2=         2~~2        '     B2=  .2.2.  'C2=         .2.2. 


s    -h 


s(s    -h   ) 


s(s    -h   ) 


-ihm(iwe.)  -hT  hm(iooe^) 

3  "  2      2'      B3  "  .2.2. 


s    -h 


s(s    -h   ) 


2       2 
(s    -h    )  (iwe   ) 

C3~  (Zh2)°-       <2> 

s(s    -h    ) 
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The  coefficients  of  the  tensor  matrix  can  be  reduced  to  the  form 


Ai  = 


(-j) 


<-A> 


g(n1)  +  inx       g(n2)  +  in2  _ 


'■"n'o1- 


A2  = 


.(—  sin  I)  (-  sin  I) 

g^)  +  ^    +   g(n2)  +  i^2  j 


(-UN£0]' 


A„  = 


(—  cos  I)  (~    COS  I)       -, 

+ 


3      [g^)  +  i^        g(J22)  +  i^2 


(-WNe0h 


Bl  = 


r(f  sinl) 


(^sinl) 


g^)  +  jsi1        g(fi2)  +  ift2 


("UN  60)j 


r(-cos     I) 


B„  = 


(-i  sin2  I)  (-4  sin2  I) 


+ 


2     |_g(nQ)  +  i^Q      g^)  +  ^      g(n2)  +  i^z 


(-wN€0,« 


r(sin  I  cos  I)  (-i  sin  I  cos  I)  (-A  sin  I  cos  I)     _, 


B3  = 


g(nQ)  +  inQ)         g(^)  +  inx 


g(n2)  +  *22 


("*Ne0)' 


Cl  = 


(y    cos  I)  (-—  COS  I) 

+ 


g^)  +  iflj       g(n2)  +  i»2 


("'Vo*' 


r(sin  I  cos  I)  (-i  sin  I  cos  I)         (-i  sin  I  cos  I)    _. 


C2  = 


g(«0)  +  ^0 


g^)  +  ^ 


g(fi2)  +  ifl2 


^n'o*' 


.(-sin     I) 


2  2 

(-i  cos     I)  (-i  cos     I) 


3       [g(£2Q)  +   i^Q        g^)  +  iS21  g(fi2)  +  ift2 


^"N^1' 


(3) 
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whe 


re    B    =  -A   ,    C    =  -A   ,    C2  =  B    ,    and  where    QQ  =  u,  Q^  =  co  +  wR, 


£2     =00-00    . 
2  H 

The  complex,    frequency-dependent  parameter,    g,    in  the  equation  of 
motion  (7)  is  found  by  integrating  over  the  velocity  distribution  for 
electron  collision  frequencies,    v(u)  =  (constant)(u),    at  each  frequency, 

o       1       Z 


g  = 


,  -     SP  3/2  *     /•  ,  \3/2  -    _-l 

4tt     P  u 


M^t)       6XP  ('GU/kT)  dU  l 


■iO.    (4) 


3  j  [V(u)  +  i  n  ]     a 

o 
which  can  be  written  in  terms  of  Dingle,    Arndt,    and  Roz[1956] 

tabulated  functions,  £c/2    M>    £3/?  M: 

«-{5[!(?)^fi)-'@«S/l(T)}54  <5> 

where  T  =  v(u)/  ft,    u  =  kT/e,    and  v(u)  =  akT/e  =  au.     Thus,    the  implica- 
tion of  the  classical  magneto-ionic  theory,    g  ~  v,    a  constant,    real 
number,    is  not  in  general  a  valid  notion,  if  an  energy  dependence 
v  =  v(u)  is  admitted. 
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The  complex  collision  frequency,     g,    can  be  introduced  into  the 
analysis  by  evaluating  the  coefficients  of  the  quartic  (Appendix  I) 

follows: 

2  2 
a4  =  **  w  eo  (C3  +  ^O*'  (6) 

a3=  -2c2^aLB3,  (7) 

a2  =  J  €q2  [Ax  +  B2  +    2C3  -  a2(B2  +   C3)  -a*  (^  +   CJ  ] 

+   iu*Q  [-A1C3  -  B2C3  -  A2  +  B2  +    2cA2  (l-a2)]     (8) 
\  =  -21**^  (A2A3  +  AlB3)  .2.262aLB3(a2.l)  (9) 


a0  =  A1B2C3"2A2A3B3  +  ^2  -  A1B32  +   A2°3  +  ^0(A1C3  +  BZ°3 
+  A1B2+  A2  +   A2  -  B2)   -  iu,0a2  (  AlC3+  A^*  A*  +A2) 

+   ^0aL(-B2C3  "A1B2+B32-A2)    "  w"  <A1+B2+C3>  " 
+  co2e2a2:(2A1  +  B2+  C3)  +  uA2a2(V  2B2+  C3) 

2  Z    2    2   /    a       -d    \  2  2   4    a  2  2    4t,  3  3X1      2  2 

"   w  60aTaL  {  Al  +  B2)  "    w  60aT  Al  "  w  e0aLB2  "  ™  60(1~a   )    ' 

(10), 
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/         5     3\ 
The  functions  t,    (x)  (  p  =  — ,  —   I  can  be  evaluated  by  an  application  of 

quadrature  techniques  to  the  integrals  [Johler  and  Walters,    I960],     or 

other  methods  described  by  Dingle,    and  Roy  [1956], 

-1 


00 


e5/2(x)  =  (f  •')  f     P5/2(x2  +   pZ)       exp  (-p)  dp, 


(11) 


*3/2(x)"'(2! 


oo        .  -1 

p         (x     +   p    )       exp  (-p)  dp, 


0 


(12) 


The  quantities    P    and    Q    are  rewritten  ; 


Li         iooe 


0  -I 


2       2         A] 

L,  1006 


0    J 


A„    ., 


a    a      +  : — 

L    T       iooe 


0    J 


A3 


0  -I 


(13) 


1-a    + 


1006 


0-IL 


fi       2      r2+     Al 
1-a       -  £     + 

L.  1006 


0    -J 


A, 


T  1006 


0     J  L 


am&  -  — 
T*       iooe 


0    J 


Q  = 


1  -  a    + 


1006 


o  J 


A      ., 

,    2_ 

a    a     + 

L,   T        ltoen 


T  iooe 


0    -I 


B 


aTC  +  — 

L.  1006 


0     -I 


1  -  a    + 


1006 


0     JL 


2         2  1 

l-aT    -  t,     +  t^ 

L  1006 


o    -i 


aTC   ""  io3T 


aT&    A3 


1006 


0    J 


(14) 


The  constant  e  or  the  factor  icoe  can  be  eliminated  from  eq  (1  -  3, 
6  -  10,  13  and  14)  for  convenience.  The  factor  was  retained  so  that 
eq  (37),  would  be  correct  dimensionally,  ia.t  J  =  current,  amperes 
per  square  meter. 
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Fig.  1  -  Coordinate  systems,   illustrating  angle  of  incidence,   §., 

magnetic  azimuth,  §  ,    magnetic  dip,    Ijof  the  magnetic  field. 
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Fig.    3   -  Electron  density-altitude  N(z),    profiles  of  the  lower  ionosphere 
illustrating  approximate  daytime-noon  (Houston)  and  local 
ionosonde  blackout  conditions. 
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FREQUENCY,   KC/S 
Fig.    9  -  Amplitude,     |T|,     of  the  reflection  coefficient,    T,    as  a  function 
of  frequency  (f  =  10  kc/s  to  1  Mc/  s)  illustrating  a  comparison  of  the 
amplitude  for  classical  and  monoenergetic  theories,   for  propagation 
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Fig.    10   -  Amplitude,     |T|,    of  the  reflection  coefficient,    T,    as  a  function 
of  frequency  (f  -  10  kc/  s  to  1  Mc/  s)  illustrating  a  comparison  of  the 
amplitude  for  classical  and  monoenergetic  theories,   for  propagation 
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Fig.    12  -Amplitude,   |  T  [ ,    of  the  reflection  coefficient,     T,    as    a 

function  of  angle  of  incidence,    $.,    for  propagation  into  the 
west,    <|>      =  270,    based  on  monoenergetic  collision  frequencies. 
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Fig.  13  -  Comparison  of  the  amplitude  j  T      |   of  the  reflection 

coefficient  for  vertical  polarization  based  on  the  classical 
constant  collision  frequency  with  the  same  reflection  coefficient 
based  on  electron  collisions  proportional  to  energy  (mono- 
energetic)  as  a  function  of  angle  of  incidence,    $.,    for  propa- 
gation into  the  east  cj>    =9°     and  propagation  into  the  west, 
d>     =  270°. 
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Fig.  15  -  Imaginary  part  of  complex  index  of  refraction,   Im  r\, 
of  the  lower  ionosphere  for  upgoing  ordinary  (o)  and 
extraordinary  (e),   propagation  components,    illustrating 
the  comparison  of  values  based  on  the  classical  magneto- 
ionic  theory  with  those  based  on  the  monoenergetic  theory. 
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Fig.  17  -  Imaginary  part  of  complex  index  of  refraction,  Im  n, 
of  the  lower  ionosphere  for  upgoing  ordinary  (o)  and 
extraordinary  (e),    propagation  components,    illustrating 
the  comparison  of  values  based  on  the  classical  magneto- 
ionic  theory  with  those  based  on  the  monoenergetic  theory, 
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Fig.   19  -  Imaginary  part  of  complex  index  of  refraction,  Im  r\,    of  the  lower 
ionosphere  for  upgoing  ordinary  (o)  and  extraordinary  (e),    propa- 
gation components,    illustrating  the  comparison  of  values  based  on 
the  classical  magneto-ionic  theory  with  those  based  on  the  mono- 
energetic  theory. 
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Fig.    20.      Real  part  of  complex  index  of  refraction,    Re  t^,    of  the  lower 
ionosphere  for  upgoing  ordinary  (o)  and  extraordinary  (e), 
propagation  components,    illustrating  the  comparison  of  values 
based  on  the  classical  magneto-ionic  theory  with  those  based 
on  the  monoenergetic  theory. 
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Fig.    21.       Imaginary  part  of  complex  index  of  refraction,  I'm  tj,     of  the 
lower  ionosphere  for  upgoing  ordinary  (o)  and  extraordinary 
(e),    propagation  components,    illustrating  the  comparison  of 
values  based  on  the  classical  magneto-ionic  theory  with  those 
based  on  the  monoenergetic  theory. 
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Fig.    22.      Real  part  of  complex  index  of  refraction,    Re  r),    of  the  lower 
ionosphere  for  upgoing  ordinary  (o)  and  extraordinary  (e), 
propagation  components,    illustrating  the  comparison  of  values 
based  on  the  classical  magneto-ionic  theory  with  those  based 
on  the  monoenergetic  theory. 
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Fig.   23  -  Imaginary  part  of  complex  index  of  refraction,  Im  r\,  of  the  lower 
ionosphere  for  upgoing  ordinary  (o)  and  extraordinary  (e),    propa- 
gation components,    illustrating  the  comparison  of  values  based  on 
the  classical  magneto-ionic  theory  with  those  based  on  the  mono- 
energetic  theory. 
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Fig.    24.      Real  part  of  complex  index  of  refraction,     Re  r],    of  the  lower 
ionosphere  for  upgoing  ordinary  (o)  and  extraordinary  (e), 
propagation  components,    illustrating  the  comparison  of  values 
based  on  the  classical  magneto -ionic  theory  with  those  based 
on  the  monoenergetic  theory. 
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Fig.    25.      Imaginary  part  of  complex  index  of  refraction,    ImT],    of  the 
lower  ionosphere  for  upgoing  ordinary  (o)  and  extraordinary 
(e),    propagation  components,    illustrating  the  comparison  of 
values  based  on  the  classical  magneto-ionic  theory  with  those 
based  on  the  monoenergetic  theory. 

4  U.  S.    GOVERNMENT  PRINTING  OFFICE  :  1963  O  -  679589 


U.  S.  DEPARTMENT  OF  COMMERCE 

Luther  H.  Hodges,  Secretary 

NATIONAL  BUREAU  OF  STANDARDS 
A.  V.  Astin,  Director 

THE  NATIONAL  BUREAU  OF  STANDARDS 

The  scope  of  activities  of  the  National  Bureau  of  Standards  at  its  major  laboratories  in  Washington,  D.C.,  and 
Boulder,  Colorado,  is  suggested  in  the  following  listing  of  the  divisions  and  sections  engaged  in  technical  work. 
In  general,  each  section  carries  out  specialized  research,  development,  and  engineering]  in  the  field  indicated  by 
its  title.    A  brief  description  of  the   activities,    and  of  the  resultant    publications,  appears  on  the  inside  ofthe 

front  cover. 

WASHINGTON,  U.  C. 

Electricity.  Resistance  and  Reactance.  Electrochemistry.  Electrical  Instruments.  Magnetic  Measurements 
Dielectrics.    High  Voltage. 

Metrology.  Photometry  and  Colorimetry.  Refractometry.  Photographic  Research.  Length.  Engineering  Metrology. 
Mass  and  Scale.    Volumetry  and  Densimetry. 

Heat.  Temperature  Physics.  Heat  Measurements.  Cryogenic  Physics.  Equation  of  State.  Statistical  Physics. 
Radiation  Physics.  X-ray.  Radioactivity.  Radiation  Theory.  High  Energy  Radiation.  Radiological  Equipment. 
Nucleonic  Instrumentation.    Neutron  Physics. 

Analytical  and  Inorganic  Chemistry.  Pure  Substances.  Spectrochemistry.  Solution  Chemistry.  Standard  Refer- 
ence Materials.    Applied   Analytical  Research,    Crystal  Chemistry. 

Mechanics.  Sound.  Pressure  and  Vacuum.  Fluid  Mechanics.  Engineering  Mechanics.  Rheology.  Combustion 
Controls. 

Polymers.    Macromolecules:     Synthesis  and   Structure.  Polymer  Chemistry.    Polymer  Physics.     Polymer  Charac- 
terisation.   Polymer  Evaluation  and  Testing.    Applied  Polymer  Standards  and  Research.    Dental  Research. 
Metallurgy.    Engineering  Metallurgy.   Microscopy  and  Diffraction.  Metal  Reactions.    Metal  Physics.     Electrolysis 
and  Metal  Deposition. 

Inorganic  Solids.  Engineering  Ceramics.  Class.  Solid  State  Chemistry.  Crystal  Growth.  Physical  Properties. 
Crystallography. 

Building  Research.     Structural  Engineering.     Fire   Research.     Mechanical  Systems.     Organic  Building  Materials. 
Codes  and   Safety  Standards.    Heat  Transfer.    Inorganic  Building  Materials.   Metallic  Bunding  Materials. 
Applied    Mathematics.     Numerical  Analysis.     Computation.     Statistical  Engineering.  Mathematical  Physics.    Op- 
erations Research. 

Data  Processing  Systems.  Components  and  Techniques.  Computer  Technology.  Measurements  Automation. 
Engineering  Applications.    Systems  Analysis. 

Atomic  Physics.  Spectroscopy.  Infrared  Spectroscopy.  Far  Ultraviolet  Physics.  Solid  State  Physics.  Electron 
Physics.    Atomic  Physics.    Plasma  Spectroscopy. 

Instrumentation.  Engineering  Electronics.  Electron  Devices.  Electronic  Instrumentation.  Mechanical  Instru- 
ments.   Basic  Instrumentation. 

Physical  Chemistry.  Thermochemistry.  Surface  Chemistry.  Organic  Chemistry.  Molecular  Spectroscopy.  Ele- 
mentary Processes.    Mass  Spectrometry.    Photochemistry  and  Radiation  Chemistry. 

Office  of  Weights  and  Measures. 

BOULDER,     COLO. 

Cryogenic  Engineering  Laboratory.  Cryogenic  Equipment.  Cryogenic  Processes.  Properties  of  Materials.  Cryo- 
genic Technical  Services. 

CENTRAL  RADIO  PROPAGATION  LABORATORY 

Ionosphere  Research  and  Propagation.  Low  Frequency  and  Very  Low  Frequency  Research.  Ionosphere  Re- 
search. Prediction  Services.  Sun-Earth  Relationships.  Field  Engineering.  Radio  Warning  Services.  Vertical 
Soundings  Research. 

Radio  Propagation  Engineering.    Data    Reduction    Instrumentation.      Radio    Noise.     Tropospheric  Measurements. 
Tropospheric  Analysis.    Propagation-Terrain  Effects.    Radio-Meteorology.    Lower  Atmosphere  Physics. 
Radio  Systems.     Applied   Electromagnetic  Theory.    High   Frequency  and   Very   High  Frequency  Research.    Fre- 
quency Utilization.    Modulation  Research.    Antenna  Research.    Radiodetermination. 

Upper  Atmosphere  and  Space  Physics.  Upper  Atmosphere  and  Plasma  Physics.  High  Latitude  Ionosphere 
Physics.    Ionosphere  and  Exosphere  Scatter.    Airglow  and  Aurora.    Ionospheric  Radio  Astronomy. 

RADIO  STANDARDS  LABORATORY 

Radio  Physics.    Radio  Broadcast  Service.    Radio  and  Microwave  Materials.    Atomic  Frequency  and  Time-Interval 

Standards.    Radio  Plasma.    Millimeter-Wave  Research. 

Circuit  Standards.    High  Frequency  Electrical  Standards.    High  Frequency  Calibration  Services.    High  Frequency 

Impedance  Standards.    Microwave  Calibration  Services.    Microwave  Circuit  Standards.  Low  Frequency  Calibration 

Services. 


NBS 


